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Abstract 



CN ■ Limitations of the Quasiparticle Random Phase Approximation (QRPA) are 

> 



studied within an exactly solvable model, with a two body interaction of 



■ Fermi type. A special attention is paid to the violation of the Pauli exclusion 

OS ■ principle (PEP) in solving the QRPA equation. A comparison of the exact 



solution, obtained by the diagonalization of a schematic nuclear Hamiltonian 



q ■ and those obtained within the standard QRPA, the renormalized QRPA, the 

QRPA with pertubative treatment of the PEP and the QRPA with exact 



consideration of the PEP, is presented. The agreement quality is judged in 
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terms of the quasiparticle number operator matrix elements in the ground 
state and in the first excited states, of the (3 transition amplitudes, of the 
Ikeda sum rule and of the nuclear matrix element for the double beta decay. 
We have found that restoring the PEP, the QRPA solutions are considerably 
stabilized and a better agreement with the exact solution is obtained. 
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I. INTRODUCTION 



The Quasiparticle Random Phase Approximation (QRPA) has been found to be a pow- 
erful method for describing many-body systems. Due to its simplicity, the proton-neutron 
QRPA is the nuclear structure method which has been most frequently used to interpret some 
nuclear structure aspects of the beta (/?) and double beta ((5/3) decay for open shell systems 
PHTDI . The QRPA provides a description of excited states by including some nucleon-nucleon 
correlations in the ground state. 

The QRPA equations are derived directly from the equation of motion. In deriving the 
QRPA equations two basic approximations are adopted: (i) The operator, which determines 
the excited state, is taken as linear superposition of two creation and two annihilation 
quasiparticle operators by considering the BCS basis as reference, (ii) The commutator of 
bifermion operators is replaced by its expectation value in the BCS ground state. This is 
usually called the "quasiboson approximation" (QBA). The QBA violates the Pauli exclusion 
principle (PEP) and this affects severely the theory. The terms which are left out by the QBA 
become more and more important when the ground state correlations are increased which 
results in a collapse of the QRPA solution. The approach based on the two approximations 
mentioned above will be conventionally called "the standard QRPA approach". 

Recently, the instability of the QRPA solution, caused by the PEP violation received 
much attention from the experts in the field. In order to improve the reliability of the 
standard QRPA description of the nuclear transitions, the renormalized version of the QRPA 
(RQRPA), which take into account the PEP in an approximate way, has been formulated 



|Tl| . |1^1 and applied to the (5 and (5(5 decay problems |13Hl5|. Indeed, the RQRPA does not 
collapse within the physical range of the interaction strength parameters. However, avoiding 
the collapse in the RQRPA a price had to be paid, namely the violation of the Ikeda sum 
rule US 

There is a constant interest in studying the physical consequences of violating the PEP 
by the QRPA solutions. Some definite conclusions can be drawn by using solvable models, 
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like, for example, the extensions p^Jlg] to proton-neutron systems of Lipkin or Moszkowski 
models [|i~9| , p0f , as they simulate the realistic cases either by analytical solutions or by a 
minimal computational effort. It is worthwhile mentioning that the study of different many- 
body approximations within schematic models was always of great interest and moreover it 
is currently considered of a major importance P^7| , pT| - P^ . pi| ] . 

The improvement of the PEP obedience, within the QRPA, can be achieved in two 
ways: (i) By a mapping technique the whole theory can be formulated in a boson picture. 
Such an approach has been outlined for the proton-neutron monopole Lipkin model in Ref. 
p2|-p6|. (ii) One can remain within the fermionic space and derive the elements of the 
QRPA equation, at least, perturbatively. Usually approximations violating the PEP are 
tested for Fermi transition since for this case the schematic models are simple and solvable 



|Tj,|22|-f26]|n] . Recently a more general solvable model appeared |27|j28| that exploits the 
SO (8) symmetry to include simultaneously the Fermi and Gamow- Teller transitions. Such 
a model is devoted to an extensive treatment of the proton-neutron pairing interaction. 

The goal of this work is to discuss some limitations of the standard QRPA approach, con- 
cerning the PEP violation. We shall follow the second possibility mentioned above and intro- 
duce new extensions of the standard QRPA approach within the proton-neutron monopole 
Lipkin model and point out some implications for realistic calculations. The newly intro- 
duced approximations will be compared with the exact results, revealing, in this way, the 
limits of the approximations. Our work completes the discussion of the Refs. [[[7[21 HH 
in which the same schematic model was considered. Indeed, for the first time, the QRPA 
solution with an exact consideration of the PEP is presented. 

The paper is organized as follows. In Sec. II, we describe the solvable model and specify 
the corresponding solution. Section III describes the standard QRPA and RQRPA within 
the chosen solvable model. In addition, new extensions of the standard QRPA approach, 
which take into account the PEP in an approximate way and exactly are introduced, re- 
spectively. In Sec. IV, the results obtained within the QRPA approaches are presented and 
compared with the exact results. Finally in Sec. V, we summarize the results and draw 



some conclusions. 



II. NUCLEAR HAMILTONIAN 

We assume a model Hamiltonian which includes a single-particle term, proton-proton 
and neutron-neutron pairing and a charge-dependent two-body interaction with particle- 
hole and particle-particle channels included: 

H = H p + H n + H res , (2.1) 



where 



H T = e T ^a\ m a Tm - G T S ] T S T (r = p,n), 

m 

H res = 2 X (3~ P + - 2 k P P + , (2.2) 



with 



<7t - i W nt 

* m 

ft = ^""^ a p m fl "m) = (P ) ) 

m 

P =X/ a pm^nm! P + = (P ), (2-3) 

m 

a' (a) being the particle creation (annihilation) operator and ~ indicating the time reversed 
states a\ m = (— l) 3 ' T-m a£_ m . 

The schematic Hamiltonian, given by Eqs. (|2.1|) - (|2.3|) , reproduces well the QRPA results 
of the realistic Hamiltonian containing G-matrix elements of the Bonn-OBEP potential for 



the beta and double beta decay transitions |32|-p4|. The strength x ( K ) of the particle- hole 



(particle-particle) interaction corresponds to the well-known parameter g p h (g pp ) commonly 
used in literature []T]-0 to parameterize the realistic ph (pp) interaction. 

Performing the Bogolyubov transformation for protons (r = p) and neutrons (r = n) 
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which defines the quasiparticle representation, and neglecting the scattering terms a p a n and 
a^ctp, the model Hamiltonian acquires the form 

H F = eC + X 1 A^A + X 2 {A^A^ + AA), (2.5) 

with 

C = Va t a +Va t a A ] = \a ] a ] ] J=0 

m m 

Ai = 40[x(«X + W X) - K (^X + 

A 2 = 4fi(x + K)u p v p u n v n . (2.6) 

For the sake of simplicity we used a single level case j p = j n = j and G p = G n = G which 
implies equal energies for protons and neutrons quasiparticles: e = e p = e n = Q G/2 and 
v i = ^jNi/2Q, Ui = \Jl — Ni/2Q with i = p,n (N p and N n are number of protons and 
neutrons, respectively). Q denotes the semi- degeneracy of the considered single level. 

The model Hamiltonian in Eqs. ( j2.5|) and (|2.6j ), resembles the Hamiltonian of the Lipkin 
model flTDfl , when Ai is taken equal to zero. We note that operators {A, A\ C} are generators 
for an SU(2) algebra. Indeed their mutual commutators are: 

[A, ^] = 1 - ^, [C, ^] = 2A\ [A, C] = 2A. (2.7) 

This model Hamiltonian is expected to account qualitatively for some features of realistic 
pn-QRPA calculations. Due to these expectations, it has been used to study the standard 
QRPA, renormalized QRPA as well as the higher order QRPA approximations for the many- 
body system, Refs. [fT^ , P^l . The salient feature of this Hamiltonian is that the stability of the 
approximate solutions can be discussed in comparison with the exact solution determined 
by diagonalizing Hp in the space of states 

\n >= (A + ) n \0 >, 0<n<2a (2.8) 

Here |0 > denotes the vacuum state for the quasiparticle operators. The matrix to be 
diagonalized can be easily calculated with the result: 



< n\Hp\n > = 2 e nm n + Aj ( m n+ i — m n + 
< n — 2\Hp\n > = A2 m n , 



where 



m r 



:< 0|A n (A t ) n |0 > = 



n! (2 ft)! 



(2 fi-n)! (2 0) r 
For n > 20, the norm overlaps m ra are vanishing. 



n m„ . 



(n < 20). 



(2.9) 



(2.10) 



III. QUASIPARTICLE RANDOM PHASE APPROXIMATION 



Another way to find an excited state for the model Hamiltonian ( j2.5|) , is to solve the 
corresponding QRPA equation. In what follows we shall briefly describe the basic ideas 
underlying this method for our solvable model. 

Within the QRPA, an excited state \Q > is created by applying a phonon creation 
operator on a state \rpa > having the properties: 



\Q >= Q ] \rpa >, Q\rpa >= 0. 
The simplest form for the phonon operator, in the fermionic space, is 

Q f = XA 1 - YA, 



(3.1) 



(3.2) 



where X and Y are called forward- and backward- going free variational amplitudes and 
satisfy the QRPA equation: 



(3.3) 



(a 






— &QRPA 


<u 


\ 


















A J 








~ u , 





where 



A = < rpa\[A, [H F , A*]]\rpa >, 
B — — < rpa\[A, [H F ,A]]\rpa >, 
U = < rpa\[A, A^]\rpa > . 
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(3.4) 



It is useful to introduce the notation: 



X = U 1/2 X, Y = U 1/2 Y, (3.5) 

A = U- 1 ' 2 AU- l /\ B = U- 1 ' 2 BU- 1 / 2 . (3.6) 
Then the QRPA eigenenergy Eqrpa and the new amplitudes X and Y are given by: 

s qrpa = (a 2 -bY 2 , 

j£ _ A + Sqrpa 



z2 ■ 



(A + £ QRPA ) 2 -B 
— -B 

Y=^- — j2 . (3.7) 

'{A + EqrpaY-B 



From the definition of the QRPA ground state \rpa > (|3.1| ), it follows that elements A, B 
and IA of the QRPA equation are function of the X and Y amplitudes. Due to this fact this 
non-linear eigenvalue problem could be solved only numerically by an iteration process. The 
functional dependence of A, B, U on X and Y is specific to the approximation scheme and 
influences crucially the final results. Below, we shall discuss, separately, several approaches. 

The standard QRPA: The simplest approximation scheme to calculate of A, B and U is 
the quasiboson approximation (QBA), which assumes [A, A^] ~< >= 1) i- e - A and 

are considered to be boson operators. Here | > denotes the uncorrelated BCS ground 
state. In this case, one finds the expressions: 

A = 2e + \ x , B = 2\ 2 , U = l, (3.8) 

which determine the excited state eigenenergy and wavefunction with normalization X 2 — 
Y 2 = 1. The drawback of this approximation scheme is the collapse of the standard 
QRPA solution within the physically acceptable interval for the nucleon-nucleon interac- 
tion strength. 

The renormalized QRPA: The renormalized QRPA (RQRPA) approach avoids the col- 
lapse of the QRPA solution for physical parameters of the nuclear Hamiltonian. Within the 
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RQRPA the commutator, [A, A'] is replaced with its expectation value in the ground state 
D =< rpa\[A, A*]\rpa > (renormalized QBA). This modifies the matrices A, B, U in the 
following way p^JT^j : 



Y 2 

A = 2tD + XxD 2 , B = 2\ 2 D 2 , U = D=(1 + —)~ 1 . (3.9) 

Note that the fermionic structure of the A, A^ operators is taken into account only in an 
approximate way. In the limit of D = 1, i.e. the \rpa > ground state is replaced by the BCS 
one | >, one gets the standard QRPA approach. 

It is worth to remark that in both the standard QRPA and the RQRPA, the elements 
A, B and U are evaluated by using some approximate schemes for the commutator [A, . 
If the commutator is exactly considered, i.e. the PEP is fulfilled, the matrices A, B and U 



take the form: 



a in \\\ i i \ ^ <rpa\C\rpa> <rpa\CC\rpa> 

A-(2e + Ax) - (6 + AO + X 1 — 2 

< rpa\A^ A\rpa > < rpa\A^A^\rpa > 
\i 2A 2 , (3.10) 



1 1 <rpa\C\rpa> < rpa\CC\rpa > 
& — _ o ) ~ A 2(/ - ttttJ p: H A 2 - 



, < rpalA^Alrpa > , < rpa\AA\rpa > 
u = 1 _ <rpa\C^pa> 



The calculation of the involved matrix elements requires the knowledge of the \rpa > ground 
state, determined by the condition in Eq. ( |3.1|) . The analytical form for \rpa > is known 
within QBA and renormalized QBA. For the phonon operator given by the Eq. ( |3.2|) , one 
obtains: 

\rpa > QBA = ne~ dAtAt \ >, d= -—, (3.13) 



2X 



where n stands for the normalization factor. 
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In general, it is not possible to find an explicit expression for \rpa > unless some addi- 
tional approximation is adopted. Fortunately, this can be achieved in the case of the solvable 
model considered in the present paper. By solving the Eq. ( |3.1| ), one obtains^ 

\rpa >^= JV£fl£)'(AUt)'| >, (3.14) 



with 



We note that if the ground state correlations are neglected in the calculation of A, B and 
IA in Eqs. (|3.10|) - (|3.12|) , i.e. the RPA vacuum \rpa > is replaced with BCS vacuum | >, one 
gets: A = (2e + Ai), B = A 2 (2 — ^) and U = 1. Obviously, the standard QRPA equation is 
recovered in the limit Q — > oo, i.e. the bifermion operators A + and A behave like bosons. 

By using the approximate \rpa > QBA ( |3.13| ) and exact \rpa > exc ( |3.15|) solutions for the 
QRPA ground state, one achieves, in fact, new extensions of the standard QRPA approach, 
namely the QRPA with the PEP (PP QRPA) included in an approximate manner and the 
QRPA with the PEP (EPP QRPA) fully fulfilled , respectively. Both methods go beyond 
the renormalized QRPA approach and require to evaluate ground state expectations values 
of the C, CC, A^A^ and A^A operators entering the expression for A, B and U as shown in 
Eqs. ( |3.10|) -( [3.12|) . This is performed free of any approximation. 

The QRPA with PEP: There is an important difference between this method and the 
RQRPA one, although both methods use the same ground state wave functions [see Eq. 
( [3. 121 )1. Indeed, within the RQRPA, the commutator of the two bifermion operators A and 
A^ is considered in approximate way while the PP QRPA takes it exactly. Indeed, the 
operators of interest have the expectation values: 



1 After the present paper was completed we learned that this equation has been derived also by 



other authors [30,31 
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< rpa\C\rpa > QBA = —4 n 2 d h 2 (d), 



QBA 'r^l^l'r 1 * QBA 

qba < rpa\CC\rpa > QBA = — 16 n 2 (d h2(d) — d 2 h^d)), 
qba < rpa\A^A^\rpa > QBA = n 2 h 2 (d), 

QBA < rpa\A^A\rpa > QBA = - n 2 ((2 - ~) d h 2 (d) + 2 £ ^(d)). (3.16) 

where the following notations have been used: 

ho (d) = <\e- dAA e- dAfA '\>=\, 

n 2 



/12(d) = < le-^AUV^I >, 

= E ^ J^l m W ~ " rfm 2 - ym 4 , (3.18) 



/ i4 (d) = <|e-^(A t )V dAtAt |>, 

= £WU + 2)U + D m ^~T m - (3 ' 19) 

We hope that this approach can be applied also for realistic calculations and within a large 
model space. Note that knowing \rpa > QBA , the QRPA matrices can be evaluated without 
the PEP violation, at least perturbatively with respect to the factor d. If the pertubative 
series is truncated to the quadratic terms in d, the resulting approach will be hereafter 
labeled by the abbreviation PP2 QRPA. 

The QRPA with exact PEP: This method can be formulated only for a solvable model 
for which the exact QRPA ground state can be analytically found. Within the EPP QRPA 
there is no violation of the PEP. However, one can not expect that the EPP QRPA solu- 
tion coincides with the exact solution for the first excited state of the nuclear Hamiltonian 
Q2.5|) . The difference is caused by the approximations incorporated in the construction of 
the operator determining the excited state [see Eq. ( ft.ll) ]. Therefore, from the direct 
comparison with the exact solution one may conclude how far the approximate description, 
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with the phonon operator of a simple structure, is from the exact picture (|3.1|) . The following 
expressions are used in elaborating the above defined procedure: 

n-i y 



exc. < rpa\AA 
e*c. < rpa\AA^ 
< rpa\C 

exc. 1 I 

exc. < rpa\CC 
exc. < rpa\A^A 



\rpa > 
\rpa > ex 
\rpa > 

I -i ex 

\rpa > ex 
\rpa > ex 



exc. exc 



. < rpa\A^A^\rpa > exc = N 2 £ W2, 

1=0 A 



JV 2 EA 2 (y) 2 Wi. 

iv 2 EA 2 (4/)(y) 2 ^, 
1=1 A 

iV 2 EA 2 (40 2 (y) 2 W, 
1=1 A 

, , exc. < rpa\C\rpa > £ 



2fi 



+ CTC . < rpa\AA^\rpa . (3.20) 



IV. RESULTS AND DISCUSSIONS 

In what follow we shall present the numerical results for the QRPA approaches described 
in the previous section and compare them with the values provided by diagonalizing Hp ( |2.5D . 
In order to continue and complete the discussion of this Hamiltonian given in Refs. fl7U2~2"f , 
we have chosen the same set of parameters as there: 

j = 9/2, Z = 4, N = 6, e = 1 MeV, (4.1) 

which determine the BCS amplitudes entering the lambdai^ parameters of the model Hamil- 
tonian Hp. Also we redefine the parameters k and x as i n Refs. [ |TT| , f2"2] ] : 

K _> K > = 2 Q K , X ^ X ' = 2 Q X - (4.2) 

The values \' = 0. and \' = 0-5 were adopted while the particle-particle strength k! was 
allowed to vary in the interval < k' < 2. Comparing the schematic calculations with the 
realistic ones, a value for k' close to unity is expected. 
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A. Excitation energies 

In Fig. |l] we plotted the dependence of the QRPA excitation energy and the first exci- 
tation energy obtained by diagonalizing Hp (bold solid line) on k', for x' — (upper figure) 
and x' = 0-5 (lower figure), respectively. Note that the standard QRPA breaks down for 
The RQRPA excitation energy remains real within the whole interval of k', although 
it deviates significantly from the exact solution beyond the breaking down point of the stan- 
dard QRPA. The effect is more evident for x' — 0.5. The PP2 QRPA and, especially, the 
PP QRPA energies reproduce quite well those of Hp, except for the values of k' approaching 
their minimum. The EPP QRPA, which take into account the PEP exactly, systematically 
overestimates the results obtained through the diagonalization of Hp. This difference might 
be attributed to the simple form of the phonon operator Q3.2| ). From this figure one 
remarks that the collapse is shifted to a large value of k' when the PEP is satisfied to a 
larger extent. For example for the PP2 QRPA, PP QRPA and EPP QRPA, the collapse 
appear at about k' = 1.80, k' = 2.50 and k' = 2.55 (%' = 0.5), respectively. This indicates 
that a real phase transition could take place in the region beyond k' = 2.5. 



B. Expectation values of the quasiparticle number operator 

In order to get additional information about the quality of different approximations, we 
calculate the expectation values of the quasiparticle number operator in the ground and first 
excited states. These are defined as follows: 

C C 
N =< rpa\ — \rpa >, AN =< rpa\Q—Q^\rpa > —N . (4-3) 

For the situations defined before, the results are: 

N = Y\ (QRPA), 
= Y 2 , (RQRPA), 
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= n 2 (-2d) h 2 (d), (PP QRPA), 

= exc. < rpa\C\rpa (EPP QRPA) (4.4) 

AA^=1 + 2F 2 , (QRPA), 
= 1 + 2F 2 , (RQRPA), 

= (X 2 + Y 2 ) + (X 2 — Y 2 — 1 — X2 + Y2 ) ® BA <rpa\C\rpa> QBA _ 
A " Y2 - QBA < rpa\CC\rpa > QBA (PP QRPA), 



An 



(X 2 + Y 2 ) + (X 2 - Y 2 - 1 - X<1 + Y ^ exo - < rpa \ C \ rpa >< 



Q ' 2 

A " V _ < rpa|CC|rpa (SPP QRPA). (4.5) 



4ft 

The expressions corresponding to the standard QRPA and RQRPA have been obtained by 
replacing the operator C/2 with its boson image, respectively. As for the remaining case we 
stay in the fermionic space and use the commutation algebra given by the Eq. (|2.7|). 



In Figs. | and |, iVo and AN, given by the above listed approximation schemes as well 
as the exact calculation are presented. One sees that the standard QRPA overestimates the 
ground state correlations near the collapse point, which reflects itself in a sudden increase of 
the average quasiparticle number. The RQRPA does not collapse at all and overestimates 
the exact result both for Nq and AiV. A distinct situation is produced by the QRPA 
approaches with PEP. There occurs an underestimation of the ground state correlations and 
the values of Nq are smaller than those corresponding to the exact solution. Practically, 
there is no difference between the PP2 QRPA and the PP QRPA up to the point where the 
former one collapses. The EPP QRPA, with an exact treatment of the PEP, provides the 
best agreement with the exact results obtained by diagonalizing Hp. We notice that none 
of the considered QRPA methods is able to reproduce the exact result for AiV which, after 
a certain point in the region k' = 1. — 1.1, is falling down, contrary to the behavior of other 
curves ||22|| . This is an indication that a more complex form for the phonon operator might 



be necessary. 
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C. Fermi (3 transition amplitudes 



We turn our attention now to transitions induced by the Fermi f3 ± operators. Neglecting 
the scattering term, as we did for the nuclear Hamiltonian Hp (see Sec. II), the Fermi (3^ 
operators, in the quasiparticle basis, take the form 



P~ = V2Q(u p v n A^ + v p u n A), (3+ = (/T)t. (4.6) 

The matrix elements of /3~ operator between ground and first excited states corresponding 
to different versions of the QRPA, are given as follows: 



<0j\/3~\rpa> = V2tt(Xu p v n + Yv p u n ) (QRPA), 



2Q(Xu p v n + Yv p u n )D 1/2 (RQRPA), 

l + d 2 (l-£)m 2 , 

2tt(Xu p v n + Yv p u n ) „ nJ (PP2 QRPA), 

1 + a z m 2 

2Qn 2 (Xu p v n + Yv p u n )(h (d) + 2^h 2 {d)) {PP QRPA), 



2Vt(Xu p v n + Yv p u n ) X 
(1 _ e,., <rpa\C\rpa> exc ^ ( EPP QRPA) ( 4 7) 

Expressions for (3 + are obtainable from the above equations by interchanging the u's and 
v's. 

In Figs. |] and |5] we examine the behavior of the (3~ and /5 + amplitudes respectively 
as functions of k' for two values for \' ( = 0., 0.5). One notes a rapid increase of /3~ and 
decrease of /3 + transition strengths for the standard QRPA and the PP2 QRPA, close to 
their collapse point. In general, the PP2 QRPA, the PP QRPA and the EPP QRPA methods 
reproduce better the trends of the exact results comparing them with the RQRPA. The best 
agreement with the exact values is achieved for the PP QRPA. Obviously these analysis 
demonstrate how important it is to have a correct treatment of the PEP, in describing the 
nuclear /3 transitions. It is worth noting the sensitivity of the /3 + amplitude to the details of 
nuclear structure wave functions. Indeed, by increasing the the particle-particle interaction 
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strength k, the matrix element of the (5 + transition operator reaches a vanishing value and 
therefore is no longer predictable. 

D. The Ikeda sum rule 

From the (3^ amplitudes one obtains, straightforwardly, the /3 strengths: 

S- = | < 0+\/3~\rpa > | 2 , S + = | < 0t\f3 + \rpa > | 2 . (4.8) 

For a ground state preserving the proton and neutron numbers in average, the two strengths 
are satisfying the Ikeda sum rule 

S- - S + = N - Z (4.9) 

where N and Z are the numbers of neutrons and protons, respectively. 

It is well known that the Ikeda Sum rule is automatically fullfiled in the standard QRPA 
and violated in the RQRPA fll6,IT||. According to the Fig. |6|, this is also true for the 



solvable model of this work. This figure shows also numerical results for the Ikeda sum 
rule predicted by the methods described so far, and compare them with the exact results. 
One notices that the exact results for Hp, do not fulfill the Ikeda sum rule and show a large 
deviation from the value of (N-Z) for k' for k' > 1. The origin of this phenomenon is expected 
to due to neglecting the scattering terms in the derivation of the assumed Hamiltonian Hp 
[|35|| . We hope that discrepancies concerning the Ikeda sum rule are substantially diminished 
by adding the contributions due to the quasiparticle operators a[a T) a\a\. + a T a T involved 
in the particle number operators. Indeed these contributions have been omitted so far, 
although their average on the ground state is not vanishing. One may conclude that the 
standard QRPA fails to reproduce the exact results for the chosen Hamiltonian. All other 
modifications of the QRPA (RQRPA, PP2 QRPA, PP QRPA, EPP QRPA) reproduce the 
trend of the exact solution, albeit the agreement with the exact results, for large value of k' 
(k' > 1) is rather poor. For non-vanishing k' the results lying closest to the exact ones are 
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those produced by PP QRPA and FPP QRPA. A reason might be the influence of higher 
excited states on the ground state induced by the diagonalization procedure. It is worthwhile 
to notice that improving the treatment of the PEP in the QRPA, a better agreement of the 
Ikeda sum rule with the (N-Z) value is achieved. 



E. Double beta decay matrix element 

In this section we shall focus our attention on the two-neutrino double beta decay mode, 
2uj3f3. Consequences of the previously presented approaches on the 2z//3/3-decay matrix 
element will be discussed. Within the solvable model considered here there is only one 
QRPA excited state and the corresponding nuclear 2z//3/3-decay matrix element takes the 
form 

M 2u = f < rpa\p-\0t >/ i < Ot\p-\rpa >j 
&QRPA + A 

Here, the states with subscripts "i" and "f are describing to the initial (A,Z) and final 
(A,Z+2) nuclei, respectively. We considered A to be equal to 0.5 MeV and performed the 
calculations for the following set of parameters JT7]: 



j = 19/2 (Z = Q, N = 14) -> (Z + 2 = 8, jV-2 = 12), x' = 0, 0.5 (4.11) 

Results corresponding to the matrix element Mjf , calculated with different approxima- 
tions, are shown in Fig. [7|, as function of k'. In addition we present also the exact results 
of Hp, considering only the contribution coming from the lowest intermediate state. One 
notices that the behaviors of the QRPA and the RQRPA curves are qualitatively similar 
to those found in the realistic calculations. The transition amplitude Mp v vanishes within 
the range 1.0 < k' < 1.5 for all calculations, including the exact ones. The values of M F U 
obtained by a better consideration of the Pauli exclusion principle (PP QRPA and EPP 
QRPA) are significantly less sensitive to the particle-particle interaction strength k'. Our 
studies show that the behavior of Mp u , as function of k' , is strongly influenced by the f3 + 
strength characterizing final nucleus. 
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V. SUMMARY AND CONCLUSIONS 



We have analyzed some limitations of the QRPA formalism in an solvable proton-neutron 
monopole Lipkin model. In addition to the standard QRPA and RQRPA, we introduced new 
extensions of the standard QRPA approach by improving (PP QRPA) and exactly (EPP 
QRPA) considering the Pauli exclusion principle^ These approaches have been used to study 
the behavior of different observables as function of the particle-particle interaction strength 

Due to the collapse of its solution for k' ~ 1.0, the standard QRPA reproduces worse 
the exact results of the nuclear Hamiltonian. Our studies show that the real collapse of 
the QRPA, usually associated to a phase transition, appears for a large value of k' (~ 2.5). 
For the PP2 QRPA, which is obtained from the PP QRPA by cutting the series in d at 
second order, the solution breakdown appears around k' = 1.7. This suggests that a better 
consideration of the PEP, within the QRPA, shifts the instability strength to a larger value. 

To shed more light on this problem we commented on the corrections induced by the 
ground state correlations, by plotting the average quasiparticle number versus k'. Our anal- 
ysis of the (3 transitions shows the sensitivity of the (3 + transition to correlations, included 
in the ground state, which violate the PEP. Concerning the Ikeda sum rule we have found 
that this is conserved neither by the exact solution of the quasiparticle nuclear Hamiltonian 
nor by the QRPA methods including the PEP. This discrepancy is very likely due to the 
structure of the model Hamiltonian. Also one expects, that including the contribution com- 
ing from the two quasiparticle and quasiparticle scattering operators entering the expressions 
of the particle number operators, the discrepancies are decreased. 



After the present work was completed we were aware about a certain overlap of the method 



EPP QRPA presented here and the self-consistent QRPA approach proposed in Ref. [31]. However 
there are also important differences such as i) the u and v coefficients are determined differently 
ii) the model Hamiltonians are different iii) our papers have distinct objectives. 
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The 2i//3/3-decay matrix element, calculated within the solvable model, is changing the 
sign when k' is increased, similar to what happens in realistic calculations. It was pointed 
out that Mp" calculated by the PP QRPA and EPP QRPA, is less sensitive to the details 
of nuclear structure than that predicted by the standard QRPA and RQRPA approaches. 

In discussing the results of the PP2 QRPA, PP QRPA and EPP QRPA we should keep 
in mind that the first two methods are just approximations of the last method. As expected, 
the PP QRPA reproduces the EPP QRPA solution better than the PP2 QRPA since it uses a 
complete pertubative expansion. The difference between the EPP QRPA and exact solution 
originates in the simple form of the phonon operator which ignores nonlinear terms like for 
example those proportional to the A + A + A and A + AA operators. To our knowledge, we are 
the first showing the inaccuracy coming exclusively from a simple form of the widely used 
phonon operator in the QRPA. An interesting point is the connection of the EPP QRPA and 
the PP QRPA to the exact solution. We note that there is no principle reason which would 
determine that the EPP QRPA results should reproduce the exact ones better than the PP 
QRPA. It turns out that the combined effect coming from the simple form of the phonon 
operator and the simple RPA wave functions, derived within quasiboson approximations, 
leads to a better agreement with exact solution for PP QRPA in some cases, e.g. for the 
lowest excitation energy (Fig. |l[) and f3 + transition amplitudes (Fig. ||]) and worse for other 
observables. For example the EPP QRPA reproduces the exact results better than PP 
QRPA when the expectation values of the quasiparticle number operators (Figs. § and [3]) 
and Ikeda Sum Rule (Fig. |6|) are calculated. 

The main conclusions of our analysis can be summarized as follows. 

New extensions of the standard QRPA with approximate (PP QRPA) and exact (EPP 
QRPA) consideration of the PEP were presented. These formalisms yield a better agreement 
with the exact results, obtained by diagonalizing the model Hamiltonian, than the standard 
QRPA and the renormalized QRPA approaches. 

The EPP QRPA results show that the collapse of the first excited states is far from 
the place where the standard QRPA breaks down, i.e., is achieved for larger values of the 
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particle-particle interaction strength. 

The comparison of the EPP QRPA results with those obtained with the exact eigenstates 
of Hp, points out the drawbacks coming from the simple structure of the QRPA phonon 
operator and suggests a range of applicability for this theory. Clearly, this analysis shows 
some limitations for the QRPA and RQRPA approaches. 

The results of the present paper support our hope that the PP QRPA approximations 
might work equally well in the case of realistic calculations for large model space, which 
are expected to offer more reliable results for the double beta decay transitions than the 
standard QRPA and the RQRPA ones. Indeed, the PP QRPA is based on the approximate 
QRPA ground state wave functions, derived within the QBA, which can be undoubtly found 
also in realistic models. This subject seems to be very interesting and therefore deserves 
further considerations. We intend to apply the PP2 QRPA method first for a realistic models 



with separable forces [32 
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FIGURES 

FIG. 1. Comparison of the lowest excitation energies resulting from the diagonalization of Hp 
(bold solid line - Hp diag.) with the standard QRPA (solid line - QRPA), the renormalized QRPA 
(dashed line - RQRPA), the QRPA with the Pauli exlusion principle up to second order in d (bold 
dashed line - PP2 QRPA), the QRPA with Pauli exclusion principle (bold dotted line - PP QRPA) 
and the QRPA with exact consideration of Pauli exclusion principle (bold dot-dashed line - EPP 
QRPA) values. 

FIG. 2. The expectation values of (half) the quasiparticle total number operator, C/2, in the 
ground state, versus k'. Conventions are the same as in Fig.l. 

FIG. 3. Differences between the expectation values of (half) quasiparticle number operator, 
C/2, in the first excited state and in the ground state as function of k'. Conventions are the same 
as in Fig.l. 

FIG. 4. The Fermi (5~ transition amplitudes between the ground and first excited states. Con- 
ventions are the same as in Fig.l. 

FIG. 5. The Fermi (3 + transition amplitudes between the ground and first excited states. Con- 
ventions are the same as in Fig.l. 

FIG. 6. The Fermi (3~ transition amplitudes between the ground and first excited states. Con- 
ventions are the same as in Fig.l. 

FIG. 7. The 2z//3/3-decay Fermi transition amplitude versus k' . The same notations as in Fig. 
1 are used. 
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